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I. INTRODUCTION 

It is now well established that quantum dots (QDs) 
behave as Kondo impurities at low temperatures 
Whereas many theoretical tools are available to ad- 
dress linear transport, the nonequilibrium regime is 
far less studied, although it is typically addressed in 
experiment.'* A solution of nonequilibrium transport 
through a 1-channel Kondo impurity was achieved in 
Ref. [5]; however exact results were obtained only for 
a specific point in the parameters space (Toulouse limit) . 
Another important development in this direction was the 
application of the Bcthe-ansatz and finding of many body 
scattering states^^"'^ Recently we found exact results for 
nonlinear transport close to a QCP in a double dot in 
series realizing the two impurity Kondo model (2IKM)ii 

The 2IKM consists of two impurity spins coupled to 
two channels of conduction electrons and, at the same 
time, interacting with each other through an exchange 
interaction K. Jones et. al.^ observed that a QCP at 
K = Kc separates a local singlet from a Kondo-screened 
phase, where Kc is of the order of the Kondo temperature 
Tr-. The exact critical behavior was found using con- 
formal field theorySiiS (CFT) and abelian bosonizatior^ii 
methods. Implications of the 2IKM for transport through 
double QDs were studied in Refs. ;i2,13T4T5T6T7T8]. 

The presence of a sharp quantum phase transition in 
the 2IKM became questionable soon after its discovery; 
in the mean field study in Ref. [Toj it was pointed out 
that the true QCP is restricted to the case of a special 
particle hole (P-H) symmetry assumed in Ref. @. This 
was confirmed by numerical renormalization group cal- 
culations.^*^ P-H symmetry breaking was later associated 
with two relevant potential scattering perturbations^— 
Thus, in real systems the critical behavior for K = Kc 
can be observed only above a certain crossover energy 
scale, denoted here as In order to obtain reliable 

predictions for QDs it is crucial to include the extra rel- 
evant perturbations associated with potential scattering 
in a real calculation. We achieved this task for a double 
QDi using the method developed by GanJ^ The finding 
of exact crossover results including P-H symmetry break- 



ing remains an open problem for the alternative proposed 
realization of the 2IKM by Zarand et. ali^ Compared to 
their QD system involving at least three leads, our system 
has only two leads making the nonequilibrium behavior 
more tractable. 

In this paper we generalize our previous results to a 
generic configuration ranging from scries to parallel QD 
attached to two leads; see Fig. ([l}. In this generic con- 
figuration transport from left to right occurs via differ- 
ent interfering paths. A particular feature of our results 
distinguishes the generic case from the series case: in 
the generic case the finite voltage conductance G[V] has 
the property G[V] — G[—V] ^ 0, leading to current rec- 
tification, similar to a diode. This effect results from 
interactions and is absent in a noninteracting Landauer 
description^?! An additional aim of this paper is to pro- 
vide important details on the calculation for the general 
series or parallel cases. 

The outline of the paper is as follows. In Sec. [H] the 
double QD system is presented and mapped to the 2IKM. 
In Sec. lIIII the conductance is calculated at the QCP using 
CFT methods, neglecting the effect of potentials scatter- 
ing. In Sec. IIVI we consider deviations from the QCP 
due to variations of K from Kc, and calculate the finite 
temperature crossover for the linear conductance using a 
mapping of the P-H symmetric 2IKM to the Ising model 
with a boundary magnetic field. We also apply this map- 
ping for the QD system proposed by Zarand et. al^ as 
a realization of the 2IKM. In Sec. |V] potential scattering 
is incorporated in the Hamiltonian close to the QCP, and 
in the crossover formula for the linear conductance. In 
Sec. IVII the full noncquilirium problem at finite voltage 
and temperature in the vicinity of the critical point is 
addressed. Sec. IVIII contains conclusions. We relegate 
details on the calculation of the nonlinear conductance 
using Keldysh Green functions (GFs) to the appendix. 



II. MODEL 

The physical system under consideration is shown 
schematically in Fig. ([T]). It consists of left (L) and right 
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FIG. 1: Schematic description of the physical system, con- 
sisting of two leads {L and R) and two quantum dots with 
effective spin-1/2. We use the convention that for finite fiux 
tiL — > tiic'*^ and t2R —* f2fle**'* where the tiL axe defined 
as tunneling amplitudes from lead {L,R) — > dot (i). 



{R) leads tunnel coupled to two quantum dots 1,2, with 
tunneling amplitudes iiL/j?, (i = 1,2). We assume that 
both dots are in the Kondo regime, with gate voltages ad- 
justed to give an odd number of electrons and the ti^/ji 
are sufficiently weak compared to the charging energy, 
U, so that charge fluctuations can be ignored. We write 
the effective spin-1/2 moments as Si and 82- We will be 
primarily interested in the case where iiL, t2R ^ tm, 
t2L so that the left lead is primarily coupled to dot 1 
and the right lead to dot 2 since only in this case will 
the QCP occur. Note that in the extreme case where 
tin = t2L = 0, this reduces to the series configuration 
analyzed, for example, in Ref. [H. The fluxes $l and 
are introduced in the triangular plaquettes as shown. 

In the standard fashiouj-- the conduction-electron 
channels that couple to the impurity are reduced to 
one-dimensional left moving Dirac fields ipiaix), where 
i = L,R and a =t,i are the lead and spin indices, re- 
spectively. We assume that a single mode in each lead 
couples to both impurities. Here we have linearized the 
conduction-electron dispersion around the Fermi level: 
tk = hvpk, where and k are measured relative to 
the Fermi level and Fermi wave number, respectively, x 
is a fictitious position variable conjugate to k. We set 
h = vp = I- 

We discuss the different terms which will appear in the 
model Hamiltonian, Eq. (|2.2p . An exchange interaction 



'-12 
U 



(2.1) 



between the impurity spins is generated by the interim- 
purity tunneling ti2. The impurity spins are also Kondo- 
coupled to the conduction-electron spin density at the 
origin: 



(repeated spin indices summed) 



In addition there are potential scattering (PS) terms oc 
ip^^"ipja (repeated spin indices summed). The system 



is driven out of equilibrium by a source drain voltage V. 
Thus, the Hamiltonian H is 

H = Ha + Hv + Ki2Si - §2+ Hk + Hps + H\ 

/OO 
dxij'^'°'idx^ja, 
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Hv ^ — dx^^'"{T^)l^:,^, 
^ J -00 
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with repeated lead and spin indices summed. To be com- 
plete one has to add the terms 

H' = V'ls/ ■ {Si X 5^2) + V-^'V^ ^,0(5^1 • 5^2). 

However, close to the QCP the first (second) term of 
H' has a similar effect as Hk (Hps). Therefore, up to a 
correction to the actual coupling constants, energy scales, 
and to the critical value of different parameters at the 
QCP, all of which we are not able to determine exactly, 
it is legitimate to drop H'. 

The Kondo interaction induces, via the Ruderman- 
Kittel-Kasuya-Yosida (RKKY) mechanism, an addi- 
tional contribution to the inter-impurity exchange, K = 
K12 + KpKKY, where 

KrkKY = 2{Si^i,S2=UHK^HK\Si^tS2=i) 
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fel>0,fc2<0 



tr{j(i)j(2)| 



With the parametrization of J^^^ given in Eq. (|2.6p . 
tr{j(i)j(2)} = 4j2sin2(26l)cos2|. Using Y.k ^ ^ ! de, 
where v is the density of states in the leads, and re- 
stricting the band width to |efc| < U, beyond which the 
effective spin description breaks down, one obtains a fer- 
romagnetic contribution 
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RKKY 
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We estimate the potential scattering amplitudes by 
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^li + ^21 
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(z-L,i?) 



tiLtipe"^^ + t2Lt2Re-''^'' 
U 



+ c jj-2 , (2.4) 

where c' is a constant factor of order 1. 

Until Sec. IVI Al we will assume the parity symmetry 



Si <-> S2, 



R. 



(2.5) 



However our results are not restricted to this will 
be discussed in Sec. (|VI Al) . Parity implies tiL = t2R = 
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^i; t2L = tin = t2. For finite flux $ = <I>i + tlie 
parity symmetry is preserved for <i>i — <i>j^. Calculating 
the Kondo couplings to second order in the tunneling 
amplitudes, under this symmetry, gives oc 

parameterize the hermitian exchange matrices by 

J(i) = J, J(2) = rVT^ J = J(l + cos(26I)t^ 
+ sin(20)[cos($/2)T"' - sin($/2)r^]), (2.6) 

where 



I arctan(t2/ti)|, J- 



U ' 



(2.7) 



Parity symmetry for the PS amplitudes implies 
T/^, Vl^ = (ImFjf = 0). We can estimate 



u 



^.cos- + c 



C/2 



(2.8) 



It is convenient to define even and odd channels "06,0 — 
^'^^'^ , in terms of which the parity transformation reads 

i/'e — > V'e, "00 ^ — "00- The most general form of 
Hk + Hps consistent with parity is 

Hk = Je0^'"^V'e/3-(5'l+52) 
+ Jo V'^°"^V'o/3-(5l + 52) 



Hps = yeV'^"'V'ea + K,V'^°"V'o 



(2.9) 



Indeed using Eq. (|2.6p we obtain Hk + Hps in this form 
with 

Je,o = ^ + ^ ^ +-Jr) ^ J[l ± 3i,,(2g) C0s($/2)], 

7 _ '-'L '^R^ '^L ■^r 



Ve.o = 

where 



J[cos(26l) -isin(26')sin($/2)] = iJ^le'"^-, 
Vi; + V^±{V^ + V^) 



(2.10) 



= - arctan(tan(20) sin(<i>/2)). 



(2.11) 



For finite flux has an imaginary part. To recover 
real coupling constants in Hk we remove this phase by 
a redefinition of the fields 



V'e ^ V'e = e-'^-Z^Ve, 



(2.12) 



In the ip'^ o basis, Hk has real coupling constants, Je, Jp 
and I Jm I , and it corresponds to the notation in Ref. [lO| . 



It is convenient to define ip[ = , ?/'2 = ^'^^ ■ 

Equivalently, the fields V'^ {j — 1,2) are related to the L- 



R basis by the rotation — {Aie' 
M 



M)iip'^, where 



V2 



As will be discussed in Sec. (| VI Ap . observability of the 
QCP in this system is restricted to the regime ti ^ t2, or 
equivalently small 9 (see Eq. (|2.7p ). In this limit the two 
impurity Kondo physics is especially transparent, since 
each QD is coupled essentially to one lead.— K can be 
tuned by means of ii2. For K ^ Kc the impurities are 
locked into a singlet, while for K = each impurity is 
Kondo-screened by the nearby lead. In the case of exact 
P-H symmetry, occurring for V^/ =0, those points in the 
if-parameters space are separated by a QCP at a critical 
value K = Kc^Tk^ 



III. CONDUCTANCE AT THE FIXED POINT 

In Ref. the conductance of a series double QD was 
calculated using the tunneling current operator. In this 
paper until Sec. (jVip we use the Kubo linear conductance 
formula written in terms of hulk current correlation func- 
tion. The reason for taking this different approach here is 
that it relates the conductance to correlation functions in 
certain field theories, that can be addressed using bound- 
ary conformal field theory or integrability methods. This 
allow us to express the conductance of double QDs de- 
scribed by the 2IKM in terms of correlation functions in 
the boundary Ising field theory. 

The linear conductance can be calculated from the 
Kubo formula 

2 pL pL poo 

G = lim hm — — ^ / dr / dr' / dre"*"^ 

L^cx,^^0 Sw(2L)2 J^^ 

x(J(r,T)J(r',0)). 

Here r is the physical coordinate; see Fig. ([2]). It should 
be distinguished from the fictitious coordinate x labeling 
the chiral fermions ipiaix). We define the chiral current 

densities in each lead jL{x) = tp"^ " {x)^La{x), jnix) — 

{x)ippa{x). The bulk current operator J(r) can be 
written as 



J(r) 



' jRi-r) r > 
- iL{-r) r<0 



It is useful to define the odd current jo{x) — jL{x) — 
Jr^x), since J_j^drJ{r) — J_j^dxjo{x)sgii{x). The con- 
ductance is given in terms of the odd current correlator, 



G = lim lim 



L pL 

dx I dy 

L J-L 



dre-"^^ 
X {joix, T)joiy, 0))sgn(a;j/)(3.1) 



The odd current joix) — 'tp^''" (r^Yjipja corresponds to 
the z— component of the fiavor current of the fermions in 
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FIG. 2: Illustration of the physical coordinate r running from 
left to right leads and the fictitious coordinate x labeling po- 
sition of the chiral fields ipia{x) in lead i = L,R. 



the L-R basis. We define the flavor current in terms of 
the fermions t/jj^ after the rotation Eq. (|2.12p . 



j-' = tp -^i'ja- (repeated indices summed) (3.2) 



The transformation Eq. (|2.12p amounts to a rotation in 
the flavor sector, 



jo = 2[cos(t>^{jfy - sin0„(/)^] 



(3.3) 



Consider the weak coupling limit J —t Q. It corre- 
sponds to a trivial boundary condition (BC) = 
0+) = Vl(x = 0-), iPr{x = 0+) = ^r{x = 0-), de- 
scribing free fermions with full reflection at the boundary. 
Also this BC makes apparent the continuity of the chiral 
fields ipL, 4'R at X = 0. Accordingly, the odd current 
correlator is given by 



{jo{x,T)jo{y,0))j=o - 



1 



TT^ (r -|- i{x — y)Y 



To calculate the odd current correlator at the nontrivial 
fixed point, we apply CFT methods and the Bose-Ising 
representation used in Ref. p^ . In this representation 
the 4 fermions V'iQ, are represented using a coset construc- 
tion in terms of three Wess-Zumino-Witten (WZW) non- 
linear a models, 5C/(2)f '"■^'^ x 5C/(2)f "''S"' x 5[/(2)f ™, 
together with a Z2 Ising model. The currents of the two 
SU{2)i a models are associated with the charge of each 
species V'iq and ?/'2a- The current of the SU{2)i model 
is associated with the total spin. 

Following Ref. [IQ] one may write down representa- 
tions of the various operators in the free fermion the- 
ory as product of charge (or isospin) bosons, the total 
spin boson, and the Ising field. The k = 2 WZW model 
has primary fields of spin j = (identity operator 1), 
j = 1/2 (fundamental field ga), and j — 1 (denoted (p). 
The k — 1 WZW model only has the identity operator 
and the j = 1/2 primary, Ha- Their scaling dimension 
is given by A = ■ The Ising model has three pri- 

mary fields: the identity operator 1 (A = 0), the Ising 
order parameter a {A — 1/16), and the energy operator 
e (A = 1/2). For example the fermion field is written in 
this representation as 





SU{2)i X 5(7(2)1 X SU{2)2 x Z2 


30(8) 




I! - 12 


i''' fipf 


(/)+ 


ihiy\h2)ie 





TABLE I; Bose-Ising- versus 50(8) Majorana- representation 
of the flavor current. 



The three factors have dimensions which add correctly 
to 1/2. The representation of other operators can be 
determined using the operator product expansion (OPE). 
For the Ising model the OPE gives 



This OPE is equivalent to that of the fc = 2 WZW model 
with the identifications <J g and e (j). 

Using the OPE, symmetry considerations, and con- 
sistency of scaling dimensions, we shall determine the 
representation of the odd current jo- The latter is 
related in Eq. (|3.3p to the flavor current operators 
{jfy and {jfy = U')*-yr ^ pij.gt consider (jf)' = 

Wia ~ V' W2a)- This is just the charge differ- 
ence between flavors, represented by /f — /|, where li 
is the SU{2y^°'^^'^^ current, [i = 1,2). For the operator 

la 



{jf)+ = tA'^^^V'L- we use Eq. (I^TD . 



'A'^ (a;)'02a(a;) oc lim 



(3.4) 



{x')go.{xyh^y^{x'yh2y{x)a{x')a{x). 

Consider the OPE of the fundamental field g x g = \ + 
(j). The operator under consideration is a spin singlet, 
ruling out (f) in the OPE. To account for the consistency 
of scaling dimensions we must have cr x cr — s- e. Hence 

^p2a o^ (/ii)^^('i2)ie- The Bose-Ising representation 
of the flavor current is summarized in the first column of 
Table m In the second column we consider an alternative 
50(8) representation introduced in Sec. IVII 

The main result of Ref. [l^ is that the nontrivial BC 
of the 2IKM at if = Kc corresponds to a change in 
the boundary condition occurring only in the Ising sec- 
tor of the theory: the nontrivial BC of the electrons 
corresponds to the free BC on the Ising chain, whereas 
the trivial BC for the electrons corresponds to the Ising 
model with a fixed boundary spiui^ We shall refer some- 
times to the BC of the full system at the nontrivial fixed 
point by "free" and at the trivial free fermion fixed point 
by "fixed" . 

The remaining sectors of the theory other than the 
Ising model remain unaffected. Correlation functions 
of factors belonging to sectors other than the Ising 
model have the form dictated by conformal invariance, 
{OAix,T)OA{y,0)) = i^^+ii^J_y))iA , where A is the scal- 
ing dimension of O. This form remains valid both at the 
trivial and nontrivial fixed points. On the other hand 
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correlation functions of fields from the Ising sector do 
depend on BC. There is a general formula for correla- 
tion function of primary operators for a BC obtained by 
fusion with a primary a^^^ 



(OA(x,r)OA(2/,0)) 



1 



xy > 



2A 



(3.5) 

Here Sj are elements of the modular S*— matrix. For the 
Ising model this is given by 

1/2 1/2 1/V2 
5*= I 1/2 1/2 -1/V2 
1/V2 -1/V2 

where the first, second and third rows and columns are 
labeled by the fields with scaling dimension 0, 1/2, 1/16, 
respectively. The change in BC in the 2IKM from trivial 
to nontrivial fixed points corresponds to fusion with the 
spin operator in the Ising sector Setting A = 1/2, a = 



1/16 we have 



1. Hence Eq. (|3.5p gives 



(e(x,T)e(y,0))/ree = (e(x, r)e(y, 0))/„ed ■ sgn(a;?;). (3.6) 

where up to a normalization factor (e(x, T)e{y, 0)) fixed oc 
T+i{x-y) ■ interpret this as a phase shift of 7r/2 

that the energy operator e{x) undergoes at x = 0. We 
proceed to evaluate the odd current correlation function. 
Since the crossed terms {{j^Y{i^)'^) vanish, we obtain 



free 

^ I 1 xy>Q 
I cos(2(/)™) xy <Q 



(3.7) 



One can use the Kubo formula Eq. ()3.1|) to calculate 
the conductance. However a calculation is unnecessary: 
Curiously, one obtains exactly the same result for the 
odd current correlation function, Eq. (|3.7p , assuming free 
fermions with partially transmitting BC, 

V'l(O^) = cos((/)m)V'L(0~) +jsin(0„)'(/'fl(O"), 



V'fl(0+) 



+\ - 



I sm 



^)V'l(O-) + cos{(j)m)^R{Q-). (3.8) 



This BC corresponds to transmission probability sin 0™ 
per spin. From Landauer formulai^ the linear conduc- 
tance at the nontrivial fixed point is 



2e2 2 
Go = — sm (/i™, 
h 



(3.9) 



where (pm is given in Eq. (|2.1ip . This gives Gq = 
for (f)m = 0, in particular for $ = 0, and for the case 
of a series QD (in this section — Q)} Also in the 
experimentally relevant regime ti 3> ^2 (see Sec. IVI Al) . 
corresponding to small [see Eq. (j2.7p ]. Go ^ 2e^//i. 

[Note however, that the actual BC at the nontriv- 
ial fixed point written in terms of the true fermions is 
very different from Eq. (|3.8[) . This is apparent from the 
vanishing of the one-particle S*— matrixiiS The auxiliary 
fermions satisfying linear boundary condition emerge in 
the S'0(8) representation that we shall use in Sec. (jVip .j 



A. T = phase diagram 

Having found the conductance at the QCP &i K = Kc, 
we shall consider the surrounding FL fixed points and 
draw a phase diagram. Here and until Sec. (jV| we con- 
. sider the P-H symmetric model. In this model charge 
transfer between the leads leading to finite current is 
allowed by the exchange interaction in Hk- The main 
role of = arctan(t2/ii) and flux $ is to modify the 
crossover scales Tk and Kc- We plot in Fig. ^ the 
phase diagram at fixed 9 as function of K and flux. 
The NFL state occurs along the curve K — Kc[^], 
where Kc ~ TkIJc, Jo, Jm] and Je, Jo, Jm depend on fiux 
through Eqs. (|2.10p . This curve is characterized by a fi- 
nite conductance G = Go at $ ^ 0, 2tt. It separates the 
K > Kc local singlet phase from the K < Kc Kondo- 
screened phase. 

The conductance vanishes in both FL phases. At 
K > Kc the system remains in its weak coupling limit, 
corresponding to weakly transmitting tunnel junctions. 
At < Kc a Kondo-screened phase is developed and 
the two channels "tpe and ipo participate in the screen- 
ing of the combined spin— 1 impurity. In the effective 
FL description both the even and odd channels acquire 
a phase shift of Sc ~ —So — tt/2. The conductance 
vanishes as a result of destructive interference between 
even and odd channels: an incoming electron from the 



left lead tp' 



V2 



scatters into the outgoing state 



V2 



--!/;£"* in the left lead, correspond- 



ing to full reflection. (The situation is reversed if a 7r/2 
phase shift occurs only in one channel). We point out 
that when we include P-H symmetry breaking the con- 
ductance is finite in the FL phases. 

In the Hamiltonian Eq. (|2.9p the condition J„i 7^ is 
required to mix the impurity singlet and triplet subspaces 
of the Hilbert space. At J,„ — the transition aX K — Kc 
corresponds to a level crossing between those subspaces. 
We point out that this special situation occurs in our 
system for the symmetric point ti = t2, {0 = 7r/4), and at 
zero flux. In this case, when K < Kc the conductance is 



G 



since the odd channel is decoupled, 5o — Q, and 



as a result of Kondo effect in the even channel 5c = it/ 2. 



IV. UNIVERSAL CROSSOVER AS FUNCTION 
OF INTER-IMPURITY INTERACTION A 

In the previous section we calculated the conductance 
at the critical value of the inter-impurity exchange in- 
teraction K = Kc and assuming P-H symmetry. In this 
situation the system flows from weak coupling (J = 0) 
to a NFL fixed point, corresponding to free BC in the 
Ising sector. At finite \K — Kc\ the system flows to an- 
other fixed points as illustrated in Fig. ([4]). Depending 
on the sign of AT — K^ those two states correspond to 
fixing the boundary spin in a semi-infinite Ising chain to 
point up or down. Note that whereas both the attrac- 
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G=0 
Kondo screening 
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FIG. 3: Schematic phase diagram as function of inter- 
impurity interaction K (controlled by ti2) and flux <& for fixed 
asymmetry (we assume a generic situation ti 7^ t2). Here P-H 
symmetry is assumed. The conductance is finite only at the 
NFL curve defined hy K = Kc[^] except for <I> = 0, 27r. 



tive and J = fixed points in Fig. ^ correspond to 
fixed boundary condition in the Ising sector, they differ 
by the impurity spin states. The latter are decoupled 
and contribute to the ground state degeneracy only at 
the repulsive J — fixed point. 



J=0 



t (fixed) 

NFL 

(free) 



K-Kc , h 



(fixed) 

FIG. 4: Schematic flow diagram of the P-H symmetric 2- 
impurity Kondo model. 



The crossover along the horizontal line in Fig. H]) is 
governed by the Kondo energy scale Tk- The universality 
of this crossover between fixed and free boundary condi- 
tions in the Ising sector, however, is spoiled by irrelevant 
operators emerging from other sectors of the theory. 

The crossover along the vertical line in Fig. (j4]) is gov- 



erned by an energy scale T* = ci 



Tk 



where Ci is 



a numerical factor of 0(1). It was argued in Ref [T3| 
that this crossover is in the universality class of the 
crossover from free- to fixed-boundary conditions in the 
Ising model, driven by a local magnetic field h at the 
boundary of a quantum Ising chain. The magnetic field 
h is linearly related to K ~ Kc, 



(K ~ 2 

T* = ci- „ = 



K 



(4.1) 



When T* <C Tk we may safely ignore sectors other than 
the Ising sector in the low energy crossover. This map- 



ping of the 2IKM to the boundary Ising model opens 
the possibility to calculate the full crossover formula for 
the conductance as function of /i ^ ^VJ^" at finite tem- 

/ J-K 



perature due to exact solvability of the boundary Ising 
model. 



A. Boundary Ising model 

It is well known that the scaling limit of the two di- 
mensional classical Ising model at its bulk critical point is 
described by a free massless Majorana field theory. Here 
we consider the two dimensional model with a boundary, 
which is equivalent to the quantum semi-infinite chain. 
After unfolding the model in the standard fashioii^^ we 
obtain a left moving Majorana fermion on the infinite 
line, 



H 



Ising 



He = hat 



(4.2) 

At /i = this model corresponds to free BC, expressed by 
the continuity of the chiral Majorana fermion field x(x) 
at a; = 0. /i is an external magnetic field acting on the 
boundary spin ctb only. Clearly h = ±00 implies fixed 
BC. The boundary spin can be written as2^ 



(4.3) 



Here a is an additional Majorana fermionic boundary de- 
gree of freedom which anticommutes with % and satisfies 
a2 = 1/2. 

The bulk energy operator of the Ising model corre- 
sponds to a mass term mxXi which is a product of a 
left- and a right-moving Majorana fields. Therefore the 
left moving factor of the energy operator, which is the 
field we refer to as the energy operator, is just the free 
Majorana fermion e{x) ~ x(^) with dimension A = 1/2. 
Note that x was introduced most naturally within free 
BC, while e was introduced to represent free fermions at 
the fixed BC fixed point. Indeed, for free BC of the Ising 
model x{x) is continuous and e{x) undergoes a 7r/2 phase 
shift at a; = [see Eq. (|3.6p ]. Hence, 



e{x) = sgn(a::)x(a;). 



(4.4) 



B. Energy correlator at finite boundary field 

In a bulk CFT a typical local operator is a product of 
left and right moving factors (t>{x) = (f)L{x)4>R{x) where 
we suppress the time variable. The Ising model has three 
primary bulk operators denoted Oa, (A = 1/2,1/16,0). 
In the presence of a boundary at a; = one can for- 
mulate the theory in terms of left moving fields only, 
(f){x) ~ 4>l{x)4>l{—x), X > Q. For example, 01/2(0;) = 
e{x)e{^x). In particular, at r = 0, y = —x the corre- 
lator of the left moving Ising fields e at any h is related 
to the one-point function of the bulk energy operator of 
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the boundary Ising model, {e{x,0)e{—x))h = {Oi/2{x))h- 
The one-point function of the bulk energy operator was 
calculated using the integrability of the boundary Ising 
model with the resuHi^^i^iiS^ 



{0„2{x))h = 



du 



Jliux 



27r 1 



(4.5) 



Here (3 — T^^ is the inverse temperature. More gen- 
erally consider the correlation function Ch{x,y,T) = 
{e{x,T)e{y,0))h- Consider a perturbative calculation of 
Ch{x,y,T) in Hb- It can be shown^^ that (i) the correc- 
tion vanishes for xy > 0, (ii) for xy < the correction 
is a function of z = t + i{x — y) . This implies that we 
can analytically continue the one-point function to find 
Ch{x,y,T), 



Ch{x,y,T) = {Oi/2{x))h\ 



x>0,y<0. (4.6) 



For X < 0,y > one can use Ch{x, y, t) = —Ch.{y, x, — r), 
where the — sign arises from the fermionic nature of e. 



Plugging this result in Eq. (|4.8p yields the result 
Giiuj^, X, y) - (zu;„, X, y) + ze""(--^) 



J2 e{sujn)0{sy)e{-sx) 



s=±l 



UJn + h^sg-n{uJn)/2 



{4.10) 



When xy > there is no dependence on h. To compare 
Eq. (|4.10|) in the nontrivial region xy < with the re- 
sult obtained by analytic continuation of the one point 
function of the energy operator, Eq. (|4.5p . we write the 
Fourier transform of Eq. (|4.10p into 



G{T,x,y) 



du e''(^+*(^-!')) ihhgn{x - y)/2 - 
ih'^sgn{x — j/)/2 



27r 



valid for x ■ y < 0. One arrives at the same result using 
Eqs. and (|i?7)) . In this notation the inte- 

gration variable u is related to the momentum of the 
particles used in the form factors method. 



C. Direct calculation of the energy correlator 



D. Finite temperature conductance 



For the present problem the desired correlator can be 
computed directly as will be done in this subsection. We 
turn to a calculation of the Majorana Green function 
(GF) Q{T,x,y) — ~{x{^i''')x{y,0)) at finite h and tem- 
perature T = . From Eq. (|4.4p . the energy correlator 
is 

{e{x, T)e{y, 0))h = -G{t, x, y)sgn{xy). (4.7) 
For h = 0, Q{T,x,y) is a free fermion GF, 

C(0)(r xy) = — 

^ ^ ' '^^ 27rsin(|(T + i(a;-y))) 



X [e{-Lo^)e{x - y) - e{Lo,,)e{y - x)], 

where w„ = ^(1 + 2n). Since the interaction in Eq. 
is quadratic in fermion fields, we may sum up the per- 
turbation series in the boundary magnetic field exactly, 
giving 

g{iujn,x,y) = g^"\iuJri,x,y) 
-f /i2^(°)(«w«, X, 0)^,(zc^„)^(°)(«w«, 0, y). (4.8) 

Here t/a(*w„) = — J^^ dre^'^'^'^ {a{T)a) is the a propaga- 
tor. When a is decoupled, its propagator is given by 
g^\iu!n) = {iujn)~^- Eq. (|4.8p becomes exact when 
Gaii'jJn) is calculated to infinite order in h. This is accom- 
plished by the self energy Ea(iwn) — ft-^tj'-^-' (zw„, 0, 0) — 
—ih?sgi\{uJn)/2. Thus 

C/a(iw„) = (zw„ -f ?ft.^sgn(w„)/2) \ (4.9) 



At finite temperature the conductance is obtained by 
analytic continuation 



G = lim lim 



le 



dx I dy sgn[xy) 



-L 



X / dre' 

1-13/2 



-(j.(x,r)jo(y,0))|_ _„„(4.11) 



where i^n 



27rn 
/3 



For h 



the finite tem- 



perature odd current correlator, {jo{x,T)jo{y,0)) free 
is given by Eq. ((3J)) where (jo(a;, r)jo(i/, 0)).7=o = 
— ■ 2r 1 / I V wT- At finite h we use Eq. (13.311 and 

(3^ sin^[^(T+i(x-y))] ^ 1 ' 

{ijO'ix,T){jf)y{y))=0, leading to 

{jo{x,T)Uy,0))H = 4cos2<^^((/)^(x,T)(j/)^(y)) 
+4sin2 (l)m{U^nx,T){jf)y{y)). 

Using the Bose Ising representation of the fiavor currents, 
given in Table HI and Eq. (|4.7p . we obtain 

{jo{x,T)joiy,0))h - -4cos2 (^W(T,a;,y))' 
-4 sin^ (f>mSgn{xy)g'^^^ (r, x, y)g (t, x, y). 

Compared to free BC, the odd current correlator obtains 
an additional term 



{io{x,T)jo{y,Q))h = {jo{x,T)jo{y,0))free 

2, G{T,x,y) - g'-°'>{T,x,y) 



-2sin^ 



/3sin(^(r + z(x-y))) 



(4.12) 



1 to the conduc- 
tance. The second term is nonvanishing only for xy < 0, 



2el 
h 
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as can be seen from Eq. (|4.10p . Note that 0{T,x,y) — 
g(o)(^r,x,y) -2^(o)(T,x,y)0(-a:y), hence 

{jo{x,T)jo{y,0))h-,oo = {jo{x,T)jo{y,0)) fixed, 

as expected. At finite T and h the contribution of the 
second term to the conductance is given by 2ReGi where 



Gi 



hm Hm - — t—ttt; / dx dy 

^{jo{x)jo{y)}[ll^^+,o+^ 



where Oo(x)jo(y))'^' = {jo{x)jo{y))h - {jo{x)jo{y)) frt 
This correlator can be expressed as a Matsubara sum 



{joix < o)jo{y > 0))i2 



4sin'^ 



/32 



/3/2 



-13/2 



rn.l 

^ Sin^ (j)m {"^T^? ^uAx-v)n 



n-1 



/3 



^ 2tt{1 + 21) + Ph^' 



The sum is evaluated as an analytic function of Vn = 
in terms of the digamma function (the logarithmic 
derivative of the gamma function) ip{z) = d\ogT{z)/dz, 



n-l 



47r 



^ 27r(l + 21) + /3h^ 



47r 2tt 



(3h^ 
47r 



Performing the analytic continuation ij/„ — > w + iO^, 
sending w — > 0, and performing the spatial integrations 
we obtain 



G/Go = l-^^[r/T*], 



F\t\ = —Re V'l 1 - + — ) , 



(4.13) 

where V'i(^) = d^\ogT{z)/dz'^ is the trigamma function 
and Go and T* are given in Eqs. and gl]). The 

scaling function F[x\ has the properties i^[0] = 1 and 
F[oo] = 0. A signature of a NFL is the existence of rele- 
vant operators in the Hamiltonian with scaling dimension 
A < 1. The QD setup discussed here allows to observe 
that the inter-impurity interaction is such a relevant per- 
turbation with A = 1/2. According to Eq. (|4.13p the 
crossover from G ^ Go to insulating FL state as function 
of if — Kc occurs at a value oi\K — Kc \ which scales with 
temperature as T^/^. 



E. Conductance in the model of Zarand et. al. 

We pause here to comment on an application of the 
Ising model with boundary magnetic field for a different 
double QD model proposed by Zarand et. alt^ as a real- 
ization of the 2IKM. We will show that in this system the 
full crossover of the conductance as function of K in the 




FIG. 5: Schematic description of the double QD proposed by 
Zarand et. al,^^ as a realization of the 2IKM. 



P-H symmetric point can be expressed in terms of the 
one point function of the spin operator of the boundary 
Ising model. 

Consider a modified QD system with an additional 
lead B coupled only to 5*2 as in Fig. Transport 
takes place between the left and right leads, where lead 
B acts as a screening channel for spin ^2. The analysis 
of Sec. [nil goes through, and the conductance is given 
by Eq. (|3.ip . where the odd current is still written as 



j^{x) = i/;t^ (T^)jV'ja, i = 1,2 = L,R. The next step 
in Sec. IIIII was to rewrite jo{x) in the basis ip' which is 
natural in the representation of the 2IKM Hamiltonian. 
For the present system this basis is 

tlL^L + tiRtpR I 

Vi = V2 = Vb- 



IR 



'*^^^/+*^/^^- is decoupled from 

V *1-L+*1H 

the impurities. We specialize to ti^ 



A third fermion ■03 



tiR. In this ba- 
sis the correlator {joix,T)jo{y,0)) occurring in the Kubo 
formula factorizes into the product of GFs for ip'i and 
-03, where the latter is a free fermion GF. For the GF of 
ip[ we use the bosonization formula Eq. (|3.4p , where the 
only factor which is sensitive to the critical point is the 
Ising spin operator, leading to 



{jo{x,T)jo{y,0)) 



(a{x,T)a{y,0)) 



TT^ (r + i(x — y)) 



2-1/ 



(4.14) 



For finite h and T Eqs. (|4.1ip and (|4.14p express the 
conductance in terms of the 2-point function for the chiral 
spin operator at finite magnetic field h. Following the 
analysis leading to Eq. (|4.6p . the 2-point function for the 
chiral spin operator is related to the one point function 
of the bulk spin operator by analytic continuation. 



(a{x,T)a{y,0))h ^ {Oi/iG{x))h\x- 



X > 0,y <0. 



The calculation of the one point function of the Ising spin 
at finite magnetic field h was addressed using integrabil- 
ity and the form factors method.— i^li^^ Different than 
the case of the energy operator, a closed expression for 
{Oi/iQ{x))h is not available. In the limiting cases h = 
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and h = ±00 CFT methods can be used.-- For BC ob- 
tained by fusion with operator a, Eq. (|3.5p gives (with 
A = 1/16) 



{a(x,T)a{y,0)) 



(T+i(2;-y))i/s y ,^ ^t-x 

Til ^ n ^ 

(r+i(2;-y))i/s -^U ^ ^ 



where Si, 



^1/16 / 5!/!^ 

'(a) = ° so/gu — • It is easy to calculate the con- 
ductance with Eq. (|4.15p . with the result G = ^ 
S'(a))- At weak coupling J — [a ~ 0) we have iSj^o) — 1; 
G = 0. This is also the result for the BC obtained by 
starting at the QCP and setting K > Kc (local singlet 



phase). At the QCP (a = 1/16) we have 



(1/16) 



0, 



G — e'^/h. In the Kondo screened phase (a — 1/2) we 
have <S'(i/2) — —1, G — 2e^ /h. We leave for a future work 
to apply Eq. (|4.14p in order to interpolate between those 
values of G at finite temperature and h cx [K — Kc). 
The additional difficulty for this system arises due to the 
presence of the a GF rather than the e GF. 



At the non-trivial fixed point the energy operator e "dis- 
appears" by double fusion; hence one obtains two rele- 
vant boundary operators (/ii)V^(/i2) and (/ii)^(/i2), with 
dimension A = 1/2. In the parity symmetric case only 
the first operator is allowed. These relevant operators 
have the dimension of a free fermion. Following Ganii 
a fermion representation emerges naturally in the SO (8) 
representation that we shall introduce in the next sub- 
section. In order for these relevant operators to have 
bosonic statistics, in the SO (8) representation indeed 
they are written as a product of a bulk fermion with 
a local fermion with dimension A = 0, which can be 
associated with a leftover impurity degree of freedom. 

On the other hand the intra-channel PS terms lead to 
marginal operators at the QCP, 



V^'IV'i ± - ii ± II 



A. Fixed point Hamiltonian in SO (8) 
representation 



(5.4) 



V. UNIVERSAL CROSSOVER AT FINITE 
POTENTIAL SCATTERING 

Until here we assumed P-H symmetry and empha- 
sized that the crossover is in the universality class of the 
boundary Ising model. Now we shall consider the more 
general situation with potential scattering (PS). We will 
see that the Ising and charge SU{2)i sectors of the the- 
ory are coupled. However this coupling can be written 
in a simple quadratic form in the Majorana S0{8) rep- 
resentation that will be introduced below. 

It is convenient to write Hps, defined in Eq. (|2.2p . in 
the V' basis [defined in Eq. ^J^], 



Hps — 



Vt 



-^(^'IV'l+V''y2) 



(5.1) 



where 



v/- -vS- p 

— - cos (j)m - ImVf; sm l 



Va = 
^ 2 



Vb - 



sin cf)^ 



' ImV^ cos ( 



(5.2) 



In the parity symmetric case Va = Vb = 0. 

At the QCP the PS terms describing charge trans- 
fer between channels V^'i and ij}' 2 generate relevant 
perturbationSfi^ To see this consider their Bose-Ising rep- 
resentation (using Table m 



'4''\'>P'2 



/i.e. - (/ii)V^(/i2)e, 
i^'\il^'^ + h.c.^ {hi)\h2)e. 



(5.3) 



Following Ref. [23| we bosonize the original theory and 
introduce four left moving bosonic fields : '■— 
^dx<j>ja- In terms of the bosons we can write the 
fermions as tjj'ja ^ FjaC"^'^^'^ . The Klein factors Fia 
take care of our sign convention required for products 
of exponentials of bosonic fields. They satisfy^ 

{F^^Fl = FlF^ - 1), {F^,F,} = 0, (5.5) 

and [-Fp, = 0, where — {«,a} and is the 
fermion number of species /i. 

Subsequently 4 linear bosonic combinations are de- 
fined, corresponding to charge, spin, fiavor, and differ- 
ence of spin between the fiavors, 

1 



2 Xl'^J"' 

ja 



ja ja 

Since the exponents of these new bosons have dimen- 
sion 1/2, we define new fermions ^ FAe~'^'^ , A = 
c,s,f,X. The new Klein factors satisfy Eq. (|5.5p with 
— c, s, /, X. To fix a convention we define^S 

pt pt — pt p p pt — pt p pt pt _ pt p 

^X^s — ^1|-'^^1J,7 ^X^s — ^2T 2J,j ^X^f ~ ^lt-'^^2t- 

The free part of the Hamiltonian can be written equiva- 
lently in bosonic or fermionic form, 

^" = E/ ^(^-'^■4)" = E/ dx^b\{^^x)^PA. 

Taking the real and imaginary parts of those fermions we 
obtain 8 Majorana fermions 



A ^ ^^A + i'A 



V2 ' 



A _ i''^A - ^A 
A2 ~ 



V2i 
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One can establish a connection between the description 

of the 2IKM in terms of 5C/(2)f '"'Sel ^ gjj^2yharge2 ^ 

SU{2)2^"^ X Z2 with 8 Majorana fermions. The two 
SU{2)i groups can be represented in terms of two bosons 

The 5C/(2)f" current f = iV/^^rfV-^s has 

the representation (j^)^ = ip\ij}s, = VSXi" V'^s- Of 

particular interest for the present work, the flavor current 
Eq. (|3.2p has the representation (see Table IJ) 

{jfy- = V^V-/, = -V2ii,^^x2- (5.6) 

The Ising fermion x can be identified with X2 ■ 
fact the nontrivial BC involves only one out of the 8 
Majorana fermions, reading X2{^ ) — ~X2^(0^)- Fo^' ^ 
description of the physics relative to the nontrivial fixed 
point it is convenient to work with the continuous Ising 
fermion field 

X{x) = sgn(x)x^ (x) = e(x)sgn(x). (5.7) 

Using Eq. (14. 3p . in the P-H symmetric case the relevant 
operator can be written as 

CTB = ix{x = 0)a = i(sgn(a;)xf (a;))x=o • a- (5.8) 

Now consider the non P-H symmetric case. From the 
SO (8) representation of the flavor current Eq. (|5.6p . the 
two PS terms in Eq. (15. 3p . (j-'^)^'^, are written in the 
trivial fixed point as 1X2 Xi 2- CFT methods tell us that 
the operators at the QCP are obtained from the operators 
at the trivial fixed point by double fusion with the spin 
operator of the Ising model. Having identified X2^ with 
the Ising fermion, double fusion gives xf^ ^ 1 + X2 ■ To 
obtain the correct bosonic statistics we argue that this 
fusion rule should be modified to 

X2 ^a + X2, 

where a is the local fermion appearing in Eq. (|5.8p . Hence 
the relevant PS operators at the QCP are 

{h,yr\h,)^ix{a, 

ih,)\h2) ^ txia. (5.9) 

Thus, a couples the Ising sector with the charge sectors. 
The main argument in favor of this form is obtained by 
considering the self-correlation function of the relevant 
operators, e.g., 

{{h%){T)(h\h,)) ^ g(°^{T,o,o)ga{T), 

at the P-H symmetric point. Fourier transforming 
Eq. (|4.9p for Gaii^^Jn) we can deduce the behavior of Ga{T): 
in the limit t <^ h'^{T ^ h'^), the correlator Ga{T) goes 
like t^{t~^). This implies that in these two limits the 
correlator ((/ij/i2)(T)(/i|/i2)) goes like t~^(t"^), respec- 
tively, as expected from an operator with scaling dimen- 
sion A = ^(1). This scaling behavior is obtained relying 



on the fact that a contains the information about the 
crossover. It explains why a, and not some other decou- 
pled local operator, should be coupled to xf and X2 
Eq. (|5.9p . On the contrary, the presence of an additional 
decoupled local operator at the QCP is ruled out as in- 
consistent with the ground state degeneracy. Away from 
the P-H symmetric point, the local operator a becomes 
also sensitive to the deviation from the QCP due to po- 
tential scattering and Ga is modified relative to Eq. (|4.9p . 

Putting together Eqs. (|5.8p and (|5.9p . the correction to 
the fixed point Hamiltonian in SO (8) representation is 

SH = i (^Aixf (a;)sgn(a;) + X2x{ix) + AsX^l^^)) 

(5.10) 

with 

K-K, 

(A2,A3) = c2VTVKReF/^,VB), (5.11) 

where and Vb are given in Eqs. (|2.8p and (|5.2p . and 
ci and C2 are constants of 0(1). This estimate of (A2, A3) 
will be justified below; as we shall see, based on the 
dimension A = 1/2 of the three relevant operators in 
Eq. (|5.10p we obtain the crossover energy scale 

T* = A^ + A^ + A^^ = A^. (5.12) 

To estimate A2 and A3 we consider the renormaliza- 
tion group flow of the inter-channel potential scattering 

I la 

operators (-0 ip'2a ^ h.c). In the presence of those 
operators the flow to the QCP stops at energy scale 
T£^. To estimate r£^ we consider the renormalization 
of these operators in the perturbative regime at energy 
scales D ^ Tk and then in the nonperturbative regime 
at energy scales D <C Tk, respectively. (A related cal- 
culation for the 2 channel Kondo model appears in |,311]). 
We assume that K — Kc. At the initial scale Dq ^ Tk 
the dimensionless bare value of these PS operators are 
Ao ~ lyRcV^ and Aq = pVb; see Eq. (|5.ip . We assume 
Ao, Ag <C 1. Since in the weak coupling regime potential 
scattering does not renormalize, we have 

AiTK) - Ao, A'{Tk) - A',. 

These can be viewed as the initial values of the coupling 
constants of the relevant perturbations (/ii)V^(/i2) and 
{hiy{h2), respectively. Since these operators have di- 
mension 1/2, the dependence of their coupling constants 
on D <^ Tk is described by 

AjP) A'{D) (TkV'^ 
A{Tk) ^ A'{Tk) J • 

The condition max{A(r£^), A' (Tps)} 1 gives the es- 
timate 

TIj, ^ ma^{TKAl Tk(A'o)2}. (5.13) 
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A more precise estimate would take into account higher 
order terms in the /9— function for A, A'. However, we 
expect that this would only change our estimate of 
by logarithmic factors. 

Identifying with \\ + in E g. ((5l^ gives the 
estimate for A2 and A3 given in Eq. (I5.11|) . Under the 
condition Aq, Aq <C 1 one has a wide energy range 
T^LR D <^Tk for the observation of the QCP. This 
can occur in a certain parameters regime, as we discuss 
in Sec. dVLa . 

We point out that our estimate for the energy scale 
T£^, Eq. (|5.13p . which agrees with [13], is inconsistent 
with that of Sakai and Shimizu, who studied the 2IKM 
with finite transfer matrix between the impurities using 
numerical renormalization group.— This discrepancy re- 
quires further investigation. 



B. Linear conductance with potential scattering 

We generalize the linear conductance calculation of 
Sec. (jlVp for finite potential scattering. Using Eqs. p.3p 
and (|5.6p the odd current operator is 



jo = 2ix2 (cos 4)mx{ + sin sgn(a;)) . (5.14) 

The operator a is now coupled to three free Majorana 
fields, and its GF, Eq. (|4.9p generalizes to Qai^'^n) = 
(iw„-|-iA^sgn(ti;„) /2) , where A^ is defined in Eq. (|5.12p . 
Similarly, 

- (Xi (a^)Xj (2/))ic^„ = (^'^n, 2;, y)Sij+hihj6g{iu;n, x, y), 
where 

dg{iujn, X, y) = Q^°'> (icj„, x, Q)ga{ii^n)Q^°'' (j^n, 0, y). 



Generalizing Eq. (|4.12p we obtain the odd current cor- 
relator 

0o(a;,r)jo(2/,0))Ai,A2,A3 = {jo{x,T)jo{y,0)) free 

+ {sm^ (/>„(A? + A^) - cos^ (j>^{Xl + Xj)) 
x4(?(°)(iw„, X, y)Sg{iuJri, x, y). 
As a result the conductance has the scaling form 

G/Go 1-F[T/T \ ^^-^ . 

(5.15) 

We see that the conductance at the free fixed point (A — 

0) is still given by Go — ^S-sin^^m- At A ^ 00 the 
Fermi liquid conductance is 

2e2 (A2)2+cos2 0„(A3)2 



G 



FL 



h A2 

We may rewrite Eq. (|5.15p as 

G — GpL 



(5.16) 



Go — Gfl 



1 - f[t/t*]. 



Gan's theory and its relation to boundary Ising 
model 



Gan presented a solution of the 2IKM, constructing 
an effective Hamiltonian for a finite region in the phase 
diagram around the critical point by controlled projec- 
tion.— The effective Hamiltoanian is solved exactly not 
only at the critical point but also for the surrounding 
Fermi-liquid phase. Excellent agreement was found with 
numerical renormalization group and CFT, in spite of 
the fact that the theory of Gan is not spin-SU(2) invari- 
ant. We shall substantiate the relation of Gan's theory 
to the CFT by showing explicitly that the operators at 
the critical point have the same form for both theories. 
In the next section we will use this approach to calculate 
the nonlinear conductance. 

Gan theory uses the 5*0(8) representation, and the 
two impurity spins turn into a local fermion d, where 



V2i 



{d,d^} = 1. Defining two Majorana fermions a 

and b — ^^j-, Gan's Hamiltonian in the P-H symmetric 
case involves only the spin-flavor (X) sector, and can be 
written as Hq = Hq^ + 5Hq where 



4"' = \ J dxx^^^.x^, 



SHg = 2iy^TKX¥ iO)b - iiK - Kc)ab. (5-17) 

We shall show that for energy scales ^ Tk this coincides 
with the Ising model Eq. (|4.2p . To see this suppose K = 
Kc and consider a mode expansion 



A 

X2{x) = ^{^pk{x)'4}k 
k 



A 

b = y^(-»fcV'fc 

k 



h.c), 



where {V'fc,V'fc'} = S{k - k'), {ipk,i^k'} — 0, and where 
initially we choose A ^ Tk as an ultraviolet cutoff. 
In the basis of tpk Gan's Hamiltonian is equal to H = 
Sfc efe^fcV'fe- One can obtain a Schrodinger's equation for 
the wave functions ipk{x) and Uk by equating the expan- 
sions of [HG,X2ix)] = [H,X2ix)] and [HG,b] = [H,b]. 



One obtains 



2i-jTK5{x)uk + idxfkix) 
-2iVTV<^fe(0) 



The solution is ipk{x) 



oc e 



ikx 



V>kiO) 



[oixH 

Uk 



Ckfkix), 
^kUk- 

(+) _L ai_ 



tan 6 



2Tk 



A-x)'fi 
- ^v^^fe(O), 

= — fc (note that we 
work with left movers). While at Tk = we have the BC 
X2"(0^) — X2"(0~), we see from the wave function that 
the effect of the first boundary term in Hg is to modify 
this BC to (0+) = -xf(0") for energies < Tk- The 
key observation is that the following operator identity 
holds if one restricts the mode expansion of its LHS and 
RHS to energies below a cutoff A <C Tk, 



=Xi(0), 



(5.18) 



K 
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where Xi{^) = X2^(2^)sgn(cc). Physically this means that 
at energy scales below Tr- the local operator b is absorbed 
into the field X2 changes its BC. Using the operator 
identity Eq. (|5.18p . we see that the term (x K — Kc in 
5Hq is equivalent to the boundary operator in the Ising 
model, Eq. (fSTS]) . 

This establishes the connection between Gan's theory 
and the boundary Ising model arising from the CFT so- 
lution, showing that Gan's anisotropic theory describes 
correctly also the vicinity of the isotropic fixed point. 

VI. CROSSOVER AT FINITE BIAS 
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FIG. 7: Scaling function for the V— asymmetric part of the 
nonlinear conductance. 



Gan's formulation of the QCP in the SO (8) Majorana 
representation provides a direct way to calculate the non- 
linear conductance at finite source drain voltage along 
the crossover from the NFL fixed point to the surround- 
ing FL fixed points, including the P-H symmetry break- 
ing. Relegating the details of the calculation based on 
the Keldysh technique to the appendix, our result is 



G 

F[t,v] 
F'[tM 



h 

Ga 

h 



Go + GsF[—, — ] + GaF'[—, —], 



1 , /I 1 

Re V"! - H h 

Ant ^ \2 Ant 2nt J 



1 

Ant 



Im ijji 



1 



1 



iv \ 

2 ' Ant ' 2ni) ' 



-\l sin^ 4'rn + A2 cos^ (/),„ + A3(l — 2 sin^ 0„j) 
A2 ■ 



sin(2(/)m)A; 



A2 sin ipm 4- Ai cos < 
' A2 



(6.1) 



Here T*,Go,(j)m and 9 are given in Eqs. (|5.12l) . dSJ]) . 
(PUj) and (EITD; iji(z) is defined below Eq. (jil^ . This 
result is valid for eV, T, T* <^Tk. When T* > T, eV the 
system is in the FL state and the nonlinear conductance 
coincides with the linear conductance, Eq. (|5.16p . Gfl = 
Go + Gs- 




FIG. 6: Scaling function for the V— symmetric part of the 
nonlinear conductance. 



The scaling functions F[t,v] and F'[t,v] are symmet- 
ric and asymmetric in v, respectively; see Figs. (jS]) and 
([7|). The asymmetric component is a new feature in the 



parallel QD, as compared to series QD where (j)m — 0. 
Having G[V] ^ G[— is a signature of interactions, 
since the Landauer noninteracting formula^ ^ leads to 
G[V] — G[—V]. This leads to a universal rectifica- 
tion effect. This rectification effect is odd under parity, 
A3 ^ — A3. Note however that it does not have a well de- 
fined transformation property with respect to $ —^ — $. 
To check the symmetry properties of our results we con- 
sidered the two impurity Anderson model for our model 
[Fig.ll])] to first order in the (intra-dot) interaction U. 
While at f/ = we have G[V] = G[-V^], which follows 
from Landauer formula, to first order in U we get a fi- 
nite G[V] — G[—V]. This asymmetric behavior of the 
conductance follows from an asymmetric dependence of 
the occupation of the dots on voltage. This simple limit 
gives the same symmetry properties of G[V] — G[—V] 
compared to the QCP, namely the rectification effect is 
odd under a parity, and does not have a well defined 
symmetry property with respect to <I> ^ — $. 

At energy scales comparable to Tk the conductance 
has additional voltage and temperature dependence due 
to irrelevant operators at the QCP. The leading irrelevant 

— 1/2 

operator is Hirr = Tj^ idxXi{^)o\x=o, with dimension 
A = 3/2.^-1° In the proposed reahzation of the 2IKM 
of Zarand et. al} '^ it leads to the conductance correc- 
tion 5G (X ^ characteristic of a NFL fixed point. 

However in the present system the irrelevant operator 
gives a nonzero correction only to fourth order, leading 

to 5G (X I as we outline below. The a-GF has an 

additional self energy = — iw^/Tj 



AT, 



1 



1 



' + iT*l2 uj + i{T*/2 + uj^/TK)' 



This has poles at oj = ~iTK^{l ± ^ 1 - For T* < 

Tk we have 



iT*/2 



K 



Qualitatively, the irrelevant correction at finite Tk has 
the same form of the fixed point conductance with 
T*/2 ^ Tk- Indeed at energy scales smaller than T*, the 
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latter has quadratic dependence on T/T* and eV jT* X- It 
should be pointed out that to fourth order in -ff^rr it is no 
longer consistent to disregard more irrelevant operators 
of dimension A = 2. However their inclusion leads only 
to the modification of the effective Kondo temperature 
in the corrections {T/TkY and {eV/TKf. 

In Figs. ([5]) and ^ the conductance is plotted in the 
parity-symmetric case at A3 = and zero temperature 
as function of source drain voltage, for different ratios 
A2/A1. The generic behavior of G\V] consists of a wide 
peak of width Tx and height Gqj with a superimposed 
narrow structure (peak or dip) of width T*, with height 
Gs (relative to the background Go). Note that Gs is pos- 
itive(negative) for Af sin^ < (>)A|cos^0m + A|(l — 
2sin^ </>m), leading to a narrow peak(dip). When A3 = 
and Aitan^m <C A2, Eq. (j6.ip predicts a peak ampli- 
tude close to the unitary limit 2e^//i. For this case, we 
mention that when and \K — Kc\ ^ , our results 
do not apply, and we expect a splitting of this peak as a 
function of V i^^^ We can obtain this behavior on a qual- 
itative level by going back to the high energy E ^ Tx 
description with Eq. (|5.17p . 



l/\/2, respectively, showing asymmetric behavior. When 
Ga > 0{Ga < 0) [defined in Eq. ([61])], the slope of the 
conductance at I^ = 0, ^\v=o, is negative(positive). 
The sign of Ga is changed under a parity transforma- 
tion (A3 — + —A3), but it also depends on the sign of the 
combination (sin(20„j)(A2 sin0„i + Ai cos </),„)). 
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FIG. 10: Nonlinear conductance at finite A3 = l/vTO for 
A2/A1 = 0, 1/4, 1/2, showing asymmetric features. 
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FIG. 8: Nonlinear conductance at A3 = for A2/A1 — 
0, 1/4, 1/2. The line shape consists of a narrow peak or dip 
structure of width T* , superimposed on top of a wide peak of 
width Tk- 




FIG. 11: Nonlinear conductance at finite A3 = l/\/2 for 
A2/A1 =0,1,5. 
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FIG. 9: Nonlinear conductance at A3 = for A2/A1 = 0, 1, 5, 
reaching the unitary limit when the relevant perturbation is 
dominated by potential scattering, namely A2 S> Ai. 



A. Observability 

In this subsection we discuss the realizability of the 
critical point in real experiment. Dealing with a repulsive 
critical point, the first condition we are concerned with 
is the smallness of the relevant perturbations, T* <C Tk- 
Secondly, we shall list some marginal corrections. 

In order to tunc K = Kc > 0, it is needed to reduce 
the ferromagnetic contribution Kukky [Eq. (|2.3p ] com- 
pared to K12 [Eq. (12. ip ]. either (i) by setting $ = tt, 
which is sufficient in the ideal situation where the device 
is perfectly parity symmetric, or (ii) in the more generic 
and realistic case, where parity symmetry is only approx- 
imate, by creating large asymmetry 



In Figs. ((To)) and (fTTj) we plot the conductance un- 
der the same conditions except A3 = 1/ \/lO and A3 = 



t 



12 



U{vJ) 



Tk± 
U vJ' 



(6.2) 
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The limit ^2 = corresponds to the the series QD. In 
either case, using Eq. (|2.ip . the condition K = ~ Tx 
is achieved by tuning the inter-dot couphng ti2 

At K = Kc (Ai 
crossover scale 



0) Eqs. ([5?T2l) and ([5?TT1) give the 



T£« = T*\k=k^ = Tk {{Rei^V^f + (lyVB)') , 

where Vb and are given in Eqs. (|5.2p and (|2.8p . 

In the parity symmetric case (i) we have Vb = 0, and 
Vjf is real and dominated by its second term in Eq. 
because $ = tt, leading to 



TIr/Tk 



U 



« 1, 



(6.3) 



as required for the validity of the critical theory. 

In the more realistic case (ii), on top of Eq. (16.2 
bound ti from below, 



we 



field, to calculate the finite temperature crossover of the 
conductance from the QCP to the stable fixed points. 
This method generalizes the CFT approach, which ad- 
dresses only the vicinity of the fixed points. We used 
this method to relate the conductance of the proposed 
system of Zarand et. al}'^ to the one-point function of 
the magnetization operator in the boundary Ising model 
which can be calculated numerically. 

Using the method developed by Gan, we solved the 
general and experimentally relevant case with potential 
scattering, and found the nonlinear conductance at finite 
temperature along the multidimensional crossover from 
QCP to surrounding FL states. Compared to the series 
double QD, we found that in the general configuration 
the universal scaling function contains both symmetric 
and asymmetric terms in the source drain voltage, lead- 
ing to a current rectification. 



VIII. ACKNOWLEDGMENTS 



such that is dominated by the first term in Eq 
In addition we demand approximate parity symmetry, 

\tlL — t2L\ 



< y/TK/Uiti/t2), sin($L - $r) < 1, 



*^^+*"-^ and t2 
1 



such that IVbI < ReV^. Here ii 

i^i^. It leads to TIJTk ^ Wlr? ~ {vJ? (|) 
1, as required for the validity of the critical theory. 

Next we estimate the marginal corrections. Spin 
SU{2) symmetry is broken by the Zeeman energy Ez = 
g^BB{Sf + S'l). This leads to a marginal operator^-'^ 
which reads in the Bose-Ising representation e. In GaAs 
QDs, the Zeeman energy is reduced due to a small g- 
factor: for the experimental conditions in Ref. 0] TV cor- 
responds to a magnetic field of few tesla, or equivalently 
to 10^ fiux quanta in a area of /xm^; for a magnetic 
field corresponding to $ = tt we have {Ez) ~ 10~^Tk, 
leading to small marginal correction to the conductance. 

Other marginal operators allowed at the QCP are the 
inter- and intra- channel PS, Eqs. (|5.3p and (|5.4p . Those 
are expected to introduce small corrections of 0((t^J)^) 
the conductance. Part of those operators break parity 
symmetry. So, the parity symmetry Eq. (|2.5p is not re- 
quired to hold exactly. Indeed the QCP has been ob- 
served numerically for a broken parity Hamiltonian in 
Ref. [13. 



VII. CONCLUSIONS 

We studied double quantum dots in the vicinity of the 
quantum critical point of the 2-impurity Kondo model. 
In the P-H symmetric model we used a mapping to 
the boundary Ising model with finite boundary magnetic 
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APPENDIX A: CALCULATION OF THE 
NONLINEAR CONDUCTANCE USING 
KELDYSH GREEN FUNCTION TECHNIQUE 

We briefiy recall basic concepts of the nonequilibrium 
formulation. Then the problem at hand will be ad- 
dressed, and the calculation of the nonlinear conductance 
will be outlined. 

One usually assumes that the system is in equilib- 
rium at some initial time, taken here to be i = —00. 
A perturbation Hi is turned on adiabatically in time, 
H = Hq + e^^Hi to drive the system out of equilibrium. 
The expectation value of an operator such as the current 
/ is given by its trace in the Heisenberg picture at t = 
weighted by the initial distribution function, 

(/) = Trie-l^^^u^O, -oo)/u(0, -00)}, 

where u(to,t) = Texp(— i J*^ dt' H{t')) and T is the 
time ordering operator. In order to employ Wick's the- 
orem, one transforms to the interaction picture, (/) = 
Tr{e~^^«u|(0, -oo)//M/(0, -cx))}, where ui{to,t) = 



aiHnt 



Oe 



iHat 



Texp{-iJ^^dt'iHi)i{t')), and Oi{t) 
Following Keldysh, for a perturbative expansion of this 
quantity it is convenient to introduce 4 types of GFs,'^- 



G"(l,l') 
Gi2(l,i') 
6-21(1,1') 
G"(l,l') 



= -z(rx(i)x(i')), 

= G<(l,l')-^(x(l')x(l)), 

= G>(l,l') = -*(x(l)x(l')), 

= -^(rx(i)x(i')>- 
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Here T is the anti-time ordering operator. It is conve- 
nient to consider an alternative set of GFs, by defining 

/ qR q< \ 

the Keldvsh GF matrix G = . , where 

- I / 



G«-^^(l,l') = T^e[±{h-m{x{l).x{l')}+). 
( Y,^ I]< \ 

Given a self energy, ^ ~ ( g 1 j the Keldysh GF 

matrix has the expansion 

G(t^) = G(°)(tj) -f G("'(^)S(w)G(°)(c^) -f (Al) 

where matrix multiplication in Keldysh space is under- 
stood and A{uj) = J dte'''^* A(t) . This leads to the Dyson 
equation for the retarded/advanced components of G, 



We see that the voltage raises the chemical potential of 
the a fermions by eV, whereas the chemical potential for 
the p fermions remains zero. The system is at equilib- 
rium at i = — oo since in this case the a and /3 fermion 
numbers are conserved. At t > —oo, Hi leads to the cur- 
rent operator / = i[Y,Hi] which drives the system out of 
equilibrium and is given by 



I — — 2i-\/TR-sin0,„Q!-)_(O)6 
-iX2 cos 0ma-|_(O)a — iX^a^ (0)a. 



(A4) 



We shall express the expectation value of the current 
by the Green functions G^^{t) where the indices refer 
to the fermion local operators = {a,b) = (1,2) and 
^1 = {a+{x = 0),a_(x = 0):P+{x = 0),/3_(x = 0)) = 
(1,2,3,4). Using Eq. (jA4p the current expectation value 
reads 



G^(cj) = G^'^"\uj) + G^^°\lu)Y^{lu)G^{lu), (A2) 
and to the Keldysh equation 

G< = G^E<G^ + (1 + G^Y^)G<^°\l + S^G^). (A3) 
We now apply this scheme to our problem with 



3 = 1 A=c,sJ,X 

eV 



Hi = SHg + «A2x{a -I- iX3X2<^' 



where Nia = J dx'ip^^°''ipia, {i — L,R). Here Hq is the 
J = fixed point Hamiltonian, including the source drain 
voltage V, and SHq is given in Eq. (|5.17p . It is more 
convenient to use SHq, which includes the local b opera- 
tor, rather than the first term in 6H, iXiX2^ {x)sgn(x)a. 
Both formulations should give the same result for energy 
scales <C Tk, as we showed generally in Sec. IV CI At 
t = —oo the system consists of two decoupled leads at 
equilibrium with different chemical potentials. It is con- 
venient to make a change of basis, in which the operator 
^ = lY^ai^La - NRa) = ^ I-ocdx^o{x) IS diagonal. 
Using Eq. (|5.14p for jo, we see that 



Y 



dxa'' a = i 



dxa- a 



where we defined new fermions a and f3: a = "'^ ^ , 
13 = ' terms of the 4 Majorana fermions ol± 

and /3± given by 

a+ = f3+ = Xi , 

a- = -(coS(/>mx{ + sin0™x?): 



-2^/T^^sin0,„G<^(^ = 0) 



(t = 0)-A3Gl_(i = 0). (A5) 



{I{t = 0)) : 

-A2 COS cj)mG'2a. 



We construct the exact GFs [appearing in Eq. (jA5p ] 
from the free GFs calculated at i = —00 {Hi = 0): for 
/i, /i' = (a+, Q!_, /?_) = (1, 2, 3, 4), one finds 



(G 



i?..A^(o) 



<^(o) 



(G 



\^ ya-|-,Q=p 



2 "mm'' 

(G<)S^^, =*/(t^), 
±]^[f{uj~eV)~ f{uj + eV)]. (A6) 



Here f{x) = (1 + e-/^)-i, j{x) = \[f{x -f eV) + f{x - 
eV)]. Note that the voltage couples the two Majorana 
fermions a+ and a_ , and here we assumed a band width 
3> V, T. The free GF for the local Majorana fermions 
V = {a,b) = (1,2) is (G«)[,"j, = 5,,'iLu + iS)-\ where S 
is a positive infinitesimal. We write Hi in a convenient 
form 

4 2 

Hi = ~iy^Xiab~i ^ ^(a+, a_, /?+, /3_)^A^^(a, 6)1., 
/I— 1 ^—1 

where 

(All, A21, A31, A41) = (-A3, A2cos(/;',„,0, -A2sin0„), 
(Ai2,A22,A32,A42) = 2 V^V (0, sin </>,„, 0, COS </>,„). (A7) 

We obtain the full GF G^^j, for z/, v' — a,b as follows. 
First suppose K = Kc (Ai = 0); we denote the different 
GFs and self energies in this case as G, S, respectively. 
At K = Kc the self energy matrix is 

S^y/ — — A^^g|°],A^/i,/, (repeated indices summed). 



6b 



13- = (sin^mXi 



cos 



<l>mX2 )■ 



Eqs. m and (|l7| give S« = -|(Ai + A^), E, 
-2»T^, E« - Sf^ = 0. Eq. m gives G« - (c. + 
^^)-\ Gf, = (c. -h 2zrK)-\ G« = Gt = 0. For 
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energies <C Tk we can approximate = (2iT/f) ^. 
For the lesser GF Eq. ((X3|) gives 



(A8) 



where matrix equation and multipUcation in ah space is 
understood. 

For K 7^ Kc the full matrix GF Gj^^/ can be calculated 
from the scries Eq. jJT]) where G(°) G and 



= 0, where acts in a6 space. Eq. (jA2p gives 
G«H ^ {u^ + ^Xy2)-\ 

G«H = (2zTK)-i-(A?/4TK)(c. + iAV2)-\ 
G^khH = T(A/2v/7V)(^ + iAV2)-^ 

For G^ , since S< = we are left with the second term of 
Eq. (|X3|l . which simphfies to [using Eq. (|X8|l and ((J2|l ] 

G< = G^S<G'^, 

where matrix equation and multiplication in ab space is 
understood. 



The GFs appearing in the current Eq. (|A5|) satisfy the 
Dyson equation 

G,^(i = 0) = 1 ^ ^ G,,,{u:)^^,,.,G^X^u). (A9) 



To evaluate G<^ we use the identity {A B)< ^ A^B"^ - 

We encounter two types of integrals for the current: 
fiu - eV) - f{LO + eV) 



I'[V,T,}?] dw 



Im ip 



1 X^/2 + ieV 



2 2ttT 
I"[V,T,X^] ^ j duj 



w + iA2/2 

, ,'1 AV2\ ^ , /I AV2 + iey 
■0 - H — - Re -0 - H ^ 



Note that = -/'[-V"], and = I"{-V\. From 

these results one can readily obtain the result for the 
nonlinear conductance, Eq. (|6.ip . 
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